Parametric amplification of vacuum fluctuations in a spinor condensate 
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Parametric amplification of vacuum fluctuations is crucial in modern quantum optics, enabling the 
creation of squeezing and entanglement. We demonstrate the parametric amplification of vacuum 
fluctuations for matter waves using a spinor F = 2 87 Rb condensate. Interatomic interactions lead 
to correlated pair creation in the rriF = ±1 states from an initial unstable tuf = condensate, 
which acts as a vacuum for vtif 0. Although this pair creation from a pure nip — condensate 
is ideally triggered by vacuum fluctuations, unavoidable spurious initial tuf = ±1 atoms induce a 
classical seed which may become the dominant triggering mechanism. We show that pair creation is 
insensitive to a classical seed for sufficiently large magnetic fields, demonstrating the dominant role 
of vacuum fluctuations. The presented system thus provides a direct path towards the generation 
of non-classical states of matter on the basis of spinor condensates. 



Parametric amplification of vacuum fluctuations in 
nonlinear media [lj plays a crucial role in the investiga- 
tion and application of non-classical states of light. These 
states have revolutionized the field of quantum optics in 
the past decades. Since the first observation of squeezed 
light [3] , these non-classical states of light have become a 
valuable tool in modern optics, e.g. for the enhancement 
of modern interferometers [31. Similarly, the production 
of entangled photon pairs [J] has triggered an ongoing 
series of fundamental tests of modern quantum mechan- 
ics [H], |(| and has many possible applications for quantum 
computing 0]. The tools developed for the production 
and manipulation of ultracold neutral atoms now bring 
many of these seminal investigations within the scope of 
experiments with matter waves. In this sense, the pro- 
duction of number-squeezed Bosc-Einstcin condensates 
(BECs) Q and spin squeezed thermal clouds has been 
recently demonstrated. 

Spinor BECs, consisting of atoms with non-zero spin 
F, provide an optimal non-linear medium for the pro- 
duction of non-classical states of matter. In these sys- 
tems, inter-particle interactions lead to a coherent popu- 
lation transfer between different Zeeman mp sublevels 
(spin dynamics). The case of a sample prepared in 
mp = (represented by |0)) is particularly interesting. 
In that case the initial stages of the spin dynamics are 
characterize d by the production of correlated atom pairs 
in |±1) (lol . resembling the production of Einstein- 
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Podolsky-Rosen (EPR) pairs in optical parametric down 
conversion 

Ideally a pure initial |0) BEC acts as a vacuum for 
atoms in mp ^ 0. Hence pair creation into |±1) can 
be understood as a parametric amplification of quantum 
vacuum fluctuations. However, parametric amplifiers are 
exponentially sensitive to any spurious initial seed in the 
amplified modes, which may easily dominate the effect of 
vacuum fluctuations. Despite careful purification procc- 
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FIG. 1: (a) Fraction of atoms transferred into \rtiF = ±1) 
within 21 ms as a function of the applied magnetic field [l2T |. 
The error bars indicate statistical uncertainties. The blue line 
is a triple Gaussian fit to guide the eye. (b) Theoretical pre- 
diction for our experimental parameters (see text), with an 
initial BEC in |0) with N = 50000 and N s = 2 seed atoms. 
The black line shows the GP result, whereas the blue line 
shows the result of our calculation, including both this clas- 
sical seed and vacuum spin fluctuations. 



dures a tiny number of spurious atoms in |±1) is exper- 
imentally unavoidable. Hence it is crucial to determine 
whether these spurious seed atoms play a significant role 
in the spin dynamics. Only in the case they do not, the 
system may safely be considered a parametric amplifier 
for vacuum fluctuations. 

In this Letter we investigate the triggering mechanism 
for the parametric amplification of correlated pairs. Sim- 
ilar to the gain of the amplification, studied in Ref. 
the triggering mechanism and its sensitivity depend cru- 
cially on the interplay between Zeeman energy, interac- 
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tions and external confinement. In particular, the rela- 
tive importance of triggering by a classical seed and by 
vacuum fluctuations depends on the magnetic field B. At 
low B the amplification is strikingly sensitive to any un- 
avoidable spurious seed. On the contrary for sufficiently 
large B the classical seed plays no significant role, and the 
system can indeed be characterized as a parametric am- 
plifier of vacuum fluctuations. The amplification of vac- 
uum fluctuations opens fascinating perspectives towards 
the analysis of two-mode squeezing and EPR entangle- 
ment of matter waves on the basis of unstable spinor 

becs [Hi EH. 

The onset of the spin dynamics is characterized by pair 
creation into |±1). In order to analyze this initial regime 
we consider the perturbed spinor i&(r, t) = (^o(r) + 
<5*(f,i))e^ t , where * (f) = (0, 0, n (r) 1/2 , 0, 0) T rep- 
resents the initial BEC in |0), u is the chemical po- 
tential, and d4?(r,t) = (<5V>-2, "^A-ii ^o, tyi, dfa) 7 ' de- 
scribes fluctuations. Up to second order in 8<p± = 
(Sijj+i ± 5$_i)/v2, pair creation is described by the 
Hamiltonian H = H+ + H- . where 



H+ = 



d 3 7 



H, 



eff + q 



(1) 

Here, q oc B 2 characterizes the quadratic Zeeman energy 
(QZE),# e// = -h 2 V 2 /2m+V trap (r)+(Ui+Uo)n (r)-n, 
Vtrap(r) is the harmonic trap, 0(f) = Uino(r), and Uq = 
(7g + 10. 92 + I8.g 4 )/35 and U x = (-7g - 5g 2 + I2.g 4 )/35 
characterize the spin-preserving and spin-changing colli- 
sions. The coupling constant associated with the s-wave 
collisional channel with total spin F is gp = 4:irh 2 ap /m. 
Note that H± is identical to the Hamiltonian describing 
an optical parametric amplifier [l3j |. 

The problem is analyzed best in the basis of eigen- 
states H e ff(p n (r) = e n ip n (r). By introducing 8<f>± = 
En^ni^n and Q ntn > = j d 3 rflip n ip n > we may rewrite 
the Hamiltonian 



h± = ( £ « + q) b. 



The Heisenberg equation ihjkdf; 



H±,Ot 



(2) 

: got 

then yields the eigenvalues ^ and the corresponding 

cigenopcrators 6± = Y<n=l (^n & «,± + R %M+n b n,±)> 
where M indicates the maximal ip n level considered. 

Imaginary eigenvalues with Im(^) > result in an ex- 
ponential amplification of spin fluctuations, which indi- 
cates the onset of pair creation. In that case, the amplifi- 
cation dynamics is dominated by the most unstable mode 
Vq with the largest imaginary part Im(£„ ) = HA(q). The 
instability rate A(g) generally shows a non-monotonous 
multi-resonant g-dependence due to the interplay be- 
tween QZE, interactions and external confinement [l2l ]. 

Based on this diagonalization of the Hamiltonian we 



can evaluate the experimentally relevant quantum evo- 
lution of the operators a n . mF and a' n mp where a n> ±i = 
(6„ j+ ± 6 ni _)/v2- The total population in |±1) is then 
P mF (t) = E„( a n,rop( < ) a i,"ii'(*)>' where the average is 
performed over the initial state \^ mF (0)) = \^>). 

In the presence of unstable spin excitation modes, 
the initial state \^) triggers the subsequent amplifica- 
tion. Ideally |^) should be a state with no particles in 
nip = ±1 (defining our vacuum state \vac}). In that 
case, pair creation is triggered by vacuum fluctuations 
(quantum triggering). However, unavoidable slight im- 
perfections always lead to N s <C N spurious atoms in 
|±f) at t = 0. Although small, these initial impurities 
serve as a classical seed of the spin dynamics (classical 
triggering), which can dominate the quantum triggering. 
It is hence crucial to determine the relative role of quan- 
tum and classical triggering. 

Note that, in our experiments [l2j the BEC in the 
|0) state is purified by quickly removing all atoms in the 
|±f) states. This results in a non-equilibrium state and 
we stress that the time after purification is too short to 
produce any thermal population of the |±1) states. Thus, 
thermal spin fluctuations do not play any role in our dis- 
cussion and impurities in |±f) can only be created by 
radio-frequency noise and magnetic field jitter after the 
purification. Since these processes act equally on each 
single atom, any spurious |±I) atom is produced in the 
same spatial mode as the original BEC. These spurious 
single-atom processes may also transfer a tiny fraction of 
the thermal cloud into |±1). However these atoms lack 
significant spatial overlap with the most unstable excita- 
tions, and hence do not contribute to the spin dynamics. 

To evaluate the triggering mechanism, the initial state 
l^) including a classical seed can be represented by 
|*) = (A^ S !)- 1 ($^ 1 <IL 1 ) JV = |vac), where ($ mjP ) ere- 
atcs (annihilates) a rap particle in the mode of the initial 
BEC. In the basis of the effective Hamiltonian, they can 
be represented by a vector $ where $ mF = J2 n ^nd n ,m F , 
with <f>„ = J d 3 TyfriQip n . Now, one obtains the total pop- 
ulation as a sum of a classically and a quantum triggered 
contribution P mF (t) = Pc(t) + -Pq(0; where 



Pcit) = AT S ** • (u\t) • U(t) + V(tf ■ V(t)j ■ (3) 
P Q {t) = Tr{v\t)-V{t)}. (4) 



the time 

,rriF (^) — 



and the operators U(t) and V(t) are 
evolution operators according to a 

En' [an,«'(*)6»',m,(0) + V 7hn ,(t)al mF (0) 

We have employed Eqs. ([3]) and ((4]) to determine the 
population at any time t. Note however that F = 2 87 Rb 
BECs present inherent hyperfine-changing losses [3], 
with an experimentally determined loss rate T ~ I0~ 2 
ms _1 . Although these losses are small (< 20% of the 
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FIG. 2: Ratio between the classically-triggered P c and 
quantum-triggered P q populations for our experimental pa- 
rameters, considering only the growth due to the most unsta- 
ble mode. The quantum triggered dynamics is characterized 
by P c < P q . 
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FIG. 3: Fraction of produced seed atoms in \m,F — ±1) as a 
function of the radio frequency power on a logarithmic scale. 
The blue line indicates a linear fit to the data which allows 
for an extrapolation to very few seed atoms. 



total number) during the typical investigation times of 
20ms, they may significantly alter the spin dynamics, 
mainly due to the dynamical modification of the reso- 
nant conditions for the instability rate A(q). We have 
carefully included these losses by splitting the evolution 
time into sub-intervals. At the beginning of each sub- 
interval we introduce losses and recalculate the evolu- 
tion operator (due to the moderate loss rate typically 
few sub-intervals already lead to convergence). Figure [T] 
(b) shows the numerical result for our experimental pa- 
rameters. Note that the multi-resonant q dependence of 
the instability rate A(q) discussed in Ref. [12j directly 
maps into a multi-peaked pair creation efficiency, which 
is in very good agreement with our previous experimen- 
tal results shown in Fig. [T] (a). As we discuss in 
detail below, the unknown average classical seed N s only 
influences the spin dynamics significantly at the low field 
resonance, and hence it is obtained from a fit to the rel- 
ative fraction of atoms on the two resonances. Our theo- 
retical results for the amplification dynamics depend on 
the total number of atoms and especially on the pre- 
cise values of the scattering lengths [l9(. We varied 
these parameters within their rather strict uncertainties 
(one standard deviation) obtaining a very good agree- 
ment with our experimental results (Fig. [I]). 

Figure Q] (b) also shows that the spin dynamics differs 
significantly from the result of a simple mean-field Gross- 
Pitaevskii (GP) approach in which field operators ipm F (^) 
are substituted by c-number fields ipm F ( , r)- Quantum 
fluctuations are absent in this description and the spin 
dynamics can only be triggered classically with ip±i(r) = 
eipo(r), with e <JC 1 (for more sophisticated approaches, 
see e.g. ref. [HI). The striking difference between GP and 
exact results indicates that the two resonances display a 
very different sensitivity to a classical seed and vacuum 
fluctuations. 

We estimate the relative importance of both triggering 
mechanisms by considering only the dominant contribu- 
tion by the most unstable mode. In that case we may 
approximate the evolution of Pc,Q{t) — Pc,Q exp(2A(g)t) 



and evaluate the time-independent ratio rj = Pc /Pq be- 
tween the classical and quantum triggering mechanisms, 
shown in Fig. [5] For sufficiently small \q\, including the 
low field resonance shown in Fig. [TJ any classical seed 
is highly relevant (77 ~ 1) due to the large overlap be- 
tween the wavefunctions of the most unstable mode and 
the original BEC. However, for larger \q\, including the 
high field resonance, the population triggered by quan- 
tum vacuum fluctuations dominates (77 -C 1) [20j. As 
detailed below, our experiments confirm this remarkable 
difference in sensitivity of the two resonances to any spu- 
rious classical seed. Note that this approach to determine 
the triggering mechanism avoids the difficulties arising 
in e.g. the investigation of fluctuations of the amplifier 
output (lij . which is largely impeded by the significant 
statistical uncertainty of the atom number measurement 
(detection noise) and small shot-to-shot variations of the 
total atom number which lead to fluctuations in the pair 
creation efficiency (amplifier noise). 

To initiate our experiments, we prepare BECs contain- 
ing 7x 10 4 F = 2 87 Rb atoms in the |0) state in an optical 
dipole trap with trapping frequencies of (176, 132, 46) Hz. 
We carefully remove residual atoms in other spin compo- 
nents by briefly applying a strong magnetic field gradient 
of w 50 G/cm, which is ramped down within 10 ms. To 
investigate the spin dynamics, the applied homogeneous 
magnetic field is subsequently lowered from 7.9 G to a 
specific value between 0.12 and 2 G within 3 ms. The 
BEC is held at the chosen magnetic field for a variable 
time to allow for spin changing collisions. Finally, the 
number of atoms in all mp components is measured in- 
dependently by applying a strong magnetic field gradient 
during time-of-fiight expansion. To evaluate the onset of 
the exponential amplification regime, we restrict our in- 
vestigation to short spin evolution times and small pop- 
ulations in |±1). 

The sensitivity of the system to a classical seed is inves- 
tigated by deliberately producing a very small symmetric 
seed population in the |±1) states prior to the spin evo- 
lution. This is accomplished by using a radio frequency 
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pulse which transfers a variable number of atoms from 
the |0) to |±1) states. This pulse is applied for 5 fis at 
a magnetic field of 7.9 G. A frequency of 5.6 MHz was 
identified, which couples the BEC in the |0) state sym- 
metrically to the |±1) states. The number of transferred 
atoms was calibrated at the smallest detectable numbers 
in the linear regime. Figure [3] shows a linear fit to the 
data which allows for an extrapolation to very small atom 
numbers. By reducing the radio-frequency power by a 
further 15 to 25 dB, it is thus possible to reproducibly 
prepare a very small number of seed atoms. Note that 
both spuriously and deliberately produced seeds result 
from similar single-atom processes, and hence have ex- 
actly the same spatial dependence as the original BEC. 
The sensitivity to the deliberately produced seed is there- 
fore representative of the sensitivity to any spurious seed 
in the experiment. 
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FIG. 4: Fraction of atoms transferred into \rtiF = ±1) as a 
function of the deliberately prepared number of seed atoms, 
(a) The fraction recorded at 0.65 G, corresponding to the 
low field resonance, shows a strong dependence on the clas- 
sical seed after an evolution time of 15 ms. (b) The fraction 
recorded after t = 23 ms on the high field resonance at 1.29 G 
is independent of the number of seed atoms and we conclude 
that it is triggered by vacuum spin fluctuations. The error 
bars indicate statistical uncertainties. The solid lines are fits 
and the dotted line is a guide to the eye indicating saturation. 



number of accidentally produced seed atoms. 

However this sensitivity to a classical seed is not gen- 
eral. Figure [4] (b) shows the fraction of transferred atoms 
for the high field resonance. For this resonance the frac- 
tion of |±1) atoms is independent of the number of de- 
liberately produced seed atoms, as expected from Fig. [21 
which indicates that a spurious classical seed is irrelevant 
for sufficiently large |g|. This experiment thus confirms 
that classical seed cannot trigger the spin dynamics in 
this mode and we conclude that pair creation indeed acts 
as a parametric amplification of vacuum fluctuations on 
the high field resonance. 

In conclusion, we have shown that a spinor F = 2 BEC 
initially prepared in an unstable |0) state can provide a 
parametric amplifier for vacuum fluctuations, where the 
|0) condensate acts as an effective vacuum for atoms in 
|±1). Similar to other parametric amplifiers the system 
is exponentially sensitive to a spurious classical seed. We 
have therefore carefully analyzed the dependence of the 
pair creation efficiency on a classical seed in |±1) . For low 
magnetic fields the amplification is extremely sensitive to 
spurious classical seed atoms, whereas for large enough 
fields the classical seed is irrelevant, and the observed pair 
creation is due to a parametric amplification of vacuum 
fluctuations. This mechanism is identical to spontaneous 
optical parametric down conversion and paves the way for 
the development of non-classical atom optics on the basis 
of unstable spinor condensates. In particular, this system 
provides a direct path towards the observation of two- 
mode squeezing of matter waves, and a promising method 
for the creation of entangled atomic EPR pairs [3, [ljj . 
In addition magnetic dipole-dipole interactions may play 
a significant role [13, IH| for the case of F = 1 BECs and 
will be the subject of future investigations. 
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Figure 2] (a) shows the fraction of atoms in the |±1) 
states after the spin dynamics depending on the number 
of seed atoms for the low field resonance (Fig. [1]). Start- 
ing at a small offset value, this fraction grows linearly 
with increasing number of seed atoms (amplification of 
23 dB) and starts to saturate for transferred fractions 
above 6%. Hence the low field resonance is strikingly 
sensitive to a classical seed, down to an extremely small 
number of seed atoms. The fact that the spin dynamics 
is sensitive to very few seed atoms also proves that the 
seed production works reliably and is not dominated by 
radio-frequency background noise. The offset is both due 
to the amplification of vacuum fluctuations and a small 
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